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Numerical simulations on the evolution of step-screw and edge-screw optical phase dislocations 
in bulk saturable self-defocusing nonlinear media are presented, with emphasis on their ability 
to induce steering waveguides for signal beams (pulses). Two schemes for directional coupling of 
such signals, both ensuring reasonable coupling-to-cross-talk efficiency ratio, are investigated. The 
parameters useful for performance optimization of the couplers are discussed. 

PACS numbers: 42.65.Tg, 42.65.Ky 
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I. INTRODUCTION 

Dark spatial solitons (DSSs) formed in bulk self- 
defocusing Kerr-type nonlinear media (NLM) are able 
to guide co-propagating probe beams/pulses pl§|| • 
The physical mechanism underlying is the intensity- 
dependent refractive-index change in a plane perpendic- 
ular to the propagation direction. Weak signal beams 
passing along these optically induced gradient waveg- 
uides H are subject to effective induced-phase modu- 
lation and are trapped. In photorefractive NLM |f| 
the situation is even more spectacular ||. Because of 
the wavelength-dependent material response, dark spa- 
tial solitons generated at low powers but at a photo- 
sensitive wavelength are able to guide and steer much 
more powerful 'signal' beams at non-photosensitive wave- 
lengths. This is confirmed for photovoltaic 0, biased 
H and quasi-steady-state |J photorefractive solitons and 
seems to hold also for incoherent dark solitons (h]] . 

Bright signal beams/pulses can be deflected by both 
electronically [|ll| and optically induced prisms fll2| . 
When they are guided by DSSs, the same can be done 
by manipulating the transverse dynamics of the dark 
beams. The transverse velocity of an optical vortex soli- 
ton (OVS) has a radial and an angular component arising 
from the transverse phase and intensity gradients, respec- 
tively Two practical ways to control the vortex 
rotation have their origin in the Guoy phase shift on both 
sides of the background beam waist [ pl|Jl5| and in the in- 
teraction of ordered structures of OVSs |16| controlled by 
the topological charges. OVS steering is demonstrated by 
superposition of a weak incoherent background field [|l7| . 
Operation of planar Y-junction splitters for signal beams 
is demonstrated in both Kerr-type jl8| and photorefrac- 
tive NLM |llj with pairs of grey DSSs born from even 
initial conditions. The possibility to branch a single in- 
put probe beam into ordered structures of sub-beams by 
quasi-two-dimensional DSSs is demonstrated numerically 



in Ref. [p0[ . Other branching and steering schemes can 
be realized by emplo yin g the inherent dynamics of r ing 
dark solitary wave s |2l| , eventual NLM saturation |22| 
and/or anisotropy p3[. 

In this work we present numerical results on how one 
input optical data channel can be linked to a selected 
output data channel by means of an optically-induced 
steering waveguide. A high steering speed of the induced 
waveguide can be obtained by using odd dark beams 
(ODBs) of finite lengths Q,|5) with a suitable choice 
of the mixed phase dislocation. Special attention is paid 
to ensure high energy efficiency of the directional coupler 
and short length of the interaction zone. Reconfigura- 
tion of the coupler is proposed to be done by changing 
the type of the phase dislocation reproduced by a multi- 
ple, active, single-voltage controlled computer-generated 
hologram p6|. 



II. STEP-SCREW AND EDGE-SCREW MIXED 
PHASE DISLOCATIONS 

The mixed phase dislocations considered consist of a 
one-dimensional phase step of limited length, which ends, 
by necessity, with pairs of phase semi-spirals with oppo- 
site helicities. In this work two possible cases are con- 
sidered. The first phase dislocation of this type, which 
will be denoted step-screw (SS), is described by the phase 
distribution 



H S p( x > tf) = A* | ^-y^sgn(y) - £ arctan 



ay 



x + b[3 



(2.1) 



where the parameters a and (3 are defined as follows: 



for \x\ < b, 

1 and /3 = — 1 for x > b, 

1 and (3=1 for x < —b. 



(2.2) 
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The phase distribution 



$ (x, y) = — — <j arctan 



x + b 



— arctan 



x — b 



(2.3) 



will be referred to as edge-screw (ES). In Eqs. ( 2. 1 )-( 2.3 ) 
the quantity A<E> stands for the magnitude of the step 
portion of the dislocation, 2b for its length, and x and y 
denote the transverse Cartesian coordinates parallel and 
perpendicular to the dislocation. Surface plots of the SS 
and ES phase dislocations are shown in Fig. and [j]b, 
respectively. 

The formation of mixed phase dislocations was first 
identified by us in interferograms of decayed crossed one- 
dimensional DSSs in the presence of moderate satura- 
tion of the nonlinearity |27[ |. From our present point of 
view the classification p4[ of such formation generated 
via an instability as an ES or SS phase dislocation can- 
not be definitive. Both distributions shown in Fig. [l] 
are characterized by the same maximal phase difference 
A$ and dislocation length 2b, but they differ in the rate 
d& / dtpl^yi, at which the phase changes azimuthally at 
their ends. In the step-screw case this change rate is twice 
as high as in the edge-screw case. In another work we 
found experimentally and confirmed numerically that in 
the interval 0.757T < A$ < 1.25w the transverse steering 
velocity of an ODB of SS type is inversely proportional 
to A$. Because of the phase gradients across the edge 
portion of the dislocation the ES one starts 'immediately' 
steering. The step portion of the SS dislocation, however, 
is forced to steer by the outlying semi-helices. The pro- 
cess stabilizes after a propagation distance of the order 
of IL^l (see Figs. 2 and 3 in Ref. |24|]), at a transverse 
velocity 2-3 times smaller than in the ES case. In view 
of this one more general statement holds: The trans- 
verse velocity of an ODB with mixed phase dislocation 
is inversely proportional to the absolute value of the az- 
imuthal phase-change rate d<& / d<p\\ x \ >b outward its ends. 
Therefore the identification of the type of ODBs gener- 
ated as a result of an instability is not very reliable. 



III. INITIAL CONDITIONS FOR THE MODEL 

The (2+l)-dimensional evolution of the ODBs in bulk 
homogeneous, isotropic and saturable NLM is described 
by the nonlinear Schrodinger equation 
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= 0, (3.1) 



where Aj_ = d 2 / 9£ 2 + d 2 / drj 2 . The transverse spatial co- 
ordinates are normalized to the initial dark beam width 
(£ = x/a, rj = y/a), and the nonlinear propagation 
path length ( is expressed in units of Rayleigh diffraction 
lengths Loiff = ka 2 . Furthermore, Lnl = (kn2lo)~ 1 
is the nonlinear length, k is the wave number inside the 



NLM, and the background beam intensity at its entrance 
isnormalized to that needed to form an ID DSS of width a 
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In order to maintain the correspondence to 



the conditions in our previous experimental work |25|,|27j , 
where NLM with saturation was used, we adopted the 
refractive-index correction 



An = n 2 \E D \ 2 /(l + s\E D 



2^7 



(3.2) 



with s = 0.3 and 7 = 3. In a weak-signal approximation 
we took into account the bright-signal beam diffraction 
and the refractive-index changes it sees as induced by the 
ODB. In the NLM its evolution is described by 
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The coupling coefficient a depends on the nature of the 
physical process of the optical nonlinearity (e.g. molec- 
ular orientation, electronic response of bound electrons, 
etc.) and on the experimental conditions (e.g. the po- 
larizations of the waves). The numerical simulations 
presented in this work refer to the wavelength ratio 
Ab/Ad ~ 1 and a — 2. The slowly-varying electric-field 
amplitudes of the ODBs were assumed to be tanh-sh&ped 
and of the form 



E D (x,y,z = 0) = V^oS(ri i0 (a;,y))tanh 



(3.4) 



where $(x,y) = $f^(x,y) for an SS, $(x,y) = 
<fr ES (x,y) for an ES phase dislocation, and r a ^(x,y) = 
\/a(x + j3b) 2 + y 2 is the effective radial coordinate. 
In order to avoid any influence of the finite back- 
ground beam of super- Gaussian form B{ri ) = 
exp{ — {\J x 1 + y 2 /w 2 ) 14 }, its width was chosen to exceed 
the maximal ODB deflection at £ = 10 more than 10 
times. The signal beam is assumed to be sec/i-shaped and 
equal in width to the ODB width a. Its intensity is chosen 
hundred times weaker than l})sg, thus not disturbin g the 
dark beams by cross-phase modulation. Equations ( jO ) 
and (3^) are solved numerically by modification of the 
beam propagation method over a 1024 x 1024 grid. 

Figure g shows the deflection of ODBs with edge-screw 
(ES) and step-screw (SS) phase dislocations vs. the non- 
linear propagation path length £ = z/Lnl for b/a = 1.0 
(solid), 2.0 (dashed), and 3.0 (dotted), normalized to the 
initial ODB width a. The difference between the trans- 
verse dynamics of the dark beams of different types is 
amazing. As already mentioned, it is due to the presence 
(in the ES case) or absence (in the SS case) of phase 
gradients across the one-dimensional portion of the dis- 
location at the entrance of the NLM. Besides by the type 
of the phase singularity, the transverse dynamics can be 
controlled in three possible ways f25| : by chaging b/a 
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ratio, the phase step A<f> and the background beam in- 
tensity. As it is seen in Fig. ||, changing the ODB length- 
to-width ratio gives an effective way to manipulate the 
ODB velocity, ft is higher for the short and slower for the 
long dislocations. The long ones however could experi- 
ence a transverse instability. More strongly pronounced 
is the increase of the velocity with decreasing the mag- 
nitude of the step portion of the phase dislocation A<E> 
p5| . but we refrained from exploiting this. The reason 
is that during the steering process both types of ODBs 
inevitably become grey (and asymptotically disappear), 
which tends to destroy the coupler modeled. The numer- 
ical calculations performed showed similar but slightly 
weaker coupling efficiencies when the transverse deflec- 
tion is controlled by A<!> than by b/a. All data in this 
work refer to A$(C = 0) = tt. Characteristic of the evo- 
lution of the ES phase dislocation is that it disturbs the 
background beam more weakly than the SS dislocation. 
In order to emphasize this, the number of gray-scale lev- 
els in the frames shown in Fig are intentionally reduced 
in the same way. The data refer to £ = z/Lpj^ = 4 and 
b/a = 1. The background beam intensity in this simula- 
tion is chosen to correspond to that of a ID DSS (Ij^sg)- 
The dashed line is intended to underline that the ES ODB 
is stronger deflected than the SS ODB. 

The background beam intensity could be additionally 
used to controll the dynamics of the ODBs, even their 
deflection does not strongly depend on it |2^] . The back- 
ground beam intensity, however, influences substantially 
the transverse beam profile and the guiding properties of 
the induced waveguide. Better spatial confinement could 
be expected at higher intensities. Numerical simulations 
were made to confirm this statement. Depending on the 
type of the dislocation (ES or SS) the highest encoupling 
efficiencies were achieved for levels (1.25 -=- 1.75)ljjg S . 
Higher intensities tend to destroy the induced waveguide 
due to the strong deformation of the dark beam profile 
under saturation of the nonlinearity. 

The main task of the numerical simulations presented 
bellow is to find the proper parameter range for b/a and 
beam intensity in order to reach maximal coupling effi- 
ciency in the analyzed schemes for directional coupling. 

IV. ONE-DIRECTIONAL COUPLER 

As a first step in this work we modeled the coupling of 
an input signal beam (pulse) to a desired output chan- 
nel by optically induced waveguides. Figure [| shows the 
positions and the forms of the weak signal beams guided 
by the corresponding dark one after a nonlinear prop- 
agation path length £ = 10. Channel is assumed to 
be supported by an OVS [^]||] with backround intensity 
1.5ljjg S , channels 1 and —1 by ODBs with SS phase dis- 
locations of b/a = 1.5 at I = 1.5ljjg S , whereas channels 
2 and —2 are considered to be induced by ODBs with ES 
phase dislocations with b/a = 2.5 at I = 1.25ljJg S . The 



change of the ODB steering direction in the NLM (i.e. 
the switching between channels with positive and nega- 
tive numbers) is controlled by reversing the gradients of 
the screw portions of the dislocations. First, real imple- 
mentation of this directional switching is feasible if there 
are reliable means to change the phase distributions of 
the ODBs. The electrically-controllable, multiply-active, 
computer-generated holograms developed j2qj are well 
suited for doing that. Because of the time needed to 
reconfigure the hologram, the relatively rare redirection 
of the data streams addressed to a desired channel will be 
preferable. Second, the energy efficiency of the guiding 
and switching should be as high as possible at the low- 
est possible cross-talk between the channels. The imag- 
inary information channels considered to be coupled to 
the exit of the NLM are assumed to be rectangle-shaped 
with a width (height) 1.5 (3.0) times the input signal- 
beam FWHM. In Fig. |] the respective positions along 
the steering direction (the arrow in Fig. ||) are sketched 
as rectangles. 

The results for the directional coupling and cross-talk 
efficiencies obtained as ratios of the transmitted energy 
in a particular channel to the total one are summarized in 
Table [jj The ODB with an SS phase profile is assumed to 
couple the input channel to channel No. 1, whereas this 
with an ES dislocation redirects the input to channel —2. 
The maximal efficiencies for channels 0, ±1, and ±2 to 
be expected are 74.7%, 70.7%, and 52.2%, respectively. 
In most channels the cross-talk signal should remain well 
below 14%. The ODB steering is accompanied by the cre- 
ation of a leading dark wave and a trailing bright peak 
on the background beam (see Fig. || in Ref. [p4| ) which 
are more strongly pronounced in the SS case. The bright 
trailing peak pushes the signal beam stronger. This qual- 
itatively explains the relatively high efficiency of the cou- 
pling to channels +1 and —1. Because of the additional 
self-defocusing of this peak it appears more effective to 
address channels +2 and —2 by ES ODBs. If the out- 
lying channels are 2.25 times larger than the input signal 
beam width a, in the ES case the efficiency in coupling 
the signal to them can be increased by 3.6% to 55.9%. 
Due to the numerical discretization the overall accuracy 
in estimating the energies is within 0.02% and the cou- 
pling efficiencies in the respective channels are correct to 
within ±0.04%. 

We should note here that because of the saturation 
of the nonlinearity the width of the induced waveguide 
is larger than in the case of pure Kerr nonlinearity. If 
a Kerr medium is considered (i.e. s = 0) the coupling 
efficiencies increase by 10%, whereas the cross-talk coef- 
ficients become slightly weaker. 

V. TWO-DIRECTIONAL COUPLER 

For addressing a desired channel from an imaginary 
array of information channels (3 x 3 in this work) we 
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used the squared-shaped channel configuration shown in 
Fig. [i| The data presented below refer to £ = 10 and 
channel width and height 1.5 times the input signal beam 
FWHM. Best spatial confinement of the signal beams is 
achieved at 1.5 times higher background beam intensity 
than that needed to form a ID DSS of width a. Be- 
sides being controlled by the type of the phase disloca- 
tion (screw, ES, SS), the ODB transverse velocity and 
deflection at the exit of the NLM are adjusted by the 
initial ODB length-to-width ratio 2b/a (b/a = 1.5 in the 
SS case; and b/a — 3.4 in the ES case; both at A<& = ir). 
Let us denote the channels as Cjj, where / and J are 
the respective row and column indices. The central out- 
put channel C22, which is on axis with respect to the 
input one, is addressed by an OVS (|||, the side chan- 
nels C21, C23, C12, and C32 by ODBs with appropriately 
±90° rotated SS phase dislocations, whereas signal beam 
coupling to each of the diagonal channels Cn, C13, C31, 
and C33 is performed by ODBs with ES phase disloca- 
tions rotated ±45°. The calculated coupling efficiencies 
are summarized in Table ||. The maximal efficiencies 
for the direct channel C22, the neighboring channel C23, 
and the diagonal channel Cn to be expected are 61.6%, 
49.6%, and 47.1% respectively, at cross-talk signals re- 
maining well below 12%. In Fig. || we composed grey- 
scale frames of the signal beams as located on the imagi- 
nary array of information channels considered. Coupling 
the input channel to information channel Cn (by an ES 
ODB) should result in undesired but equal cross-signals 
in channels Cjj and Cjj (I ^ J). If the input signal is 
coupled to channel C23 the same should hold for chan- 
nels Cu and C3J. The data obtained indicate that the 
discretization in our simulations leads to a maximal in- 
accuracy in estimating the coupling efficiencies that is 
within ±0.04%. If pure Kerr nonlinearity is considered 
the coupling efficiencies are 1% to 8% higher. 

VI. CONCLUSION 

The numerical data on the nonlinear evolution of 
mixed step-screw and edge-screw phase dislocations un- 
der equivalent conditions shows an amazing difference 
between the transverse steering dynamics of SS and ES 
ODBs. Such odd dark beams are able to induce steer- 
ing conduits for weak information beams (pulses) in bulk 
nonlinear media. Due to this two different schemes for 
directional coupling are proposed. The relatively high 
coupling and moderate cross-talk efficiencies in both of 
them provide a reasonable basis for further optimization 
and may open the way to constructing parallel rcconfig- 
urablc all-optical logical elements. 
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FIG. 1. (a) step-screw and (b) edge-screw mixed phase dis- 
locations described by Eqs. (2.1) and (2.3), respectively. 




FIG. 2. Relative deflection of ODBs with edge-screw (ES) 
and step-screw (SS) phase dislocations vs. nonlinear prop- 
agation path length £ = z/Lml for b/a = 1.0 (solid), 2.0 
(dashed), and 3.0 (dotted) at A$ = ir. 




FIG. 3. Grey-scale plots of ES and SS ODBs with initial 
length-to- width ratio b/a = 1 after nonlinear propagation 
path length f = 4. Some 3% of the total computational area 
is shown. See text for details. 
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FIG. 4. Transverse profiles (along the steering direction y) 
of the signal beams coupled to one of the five output chan- 
nels considered. The intensities are normalized to that of the 
signal in the direct channel No. 0. 
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FIG. 5. Composition of grey-scale images of the signal 
beams at the exit of the two-directional coupler. Each channel 
is addressed separately and the coupler length is £ = 10. 
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TABLE I. One-directional coupler: Efficiencies (in %) in 
addressing channels 0, 1, and —2 by an OVS, SS ODB, and 
ES ODB, respectively, and cross-talk efficiencies. 



Channel 
guiding by 


-2 


-f 





1 


2 


Energy 
losses 


OVS 


0.3 


9.5 


74.7 


9.5 


0.3 


< 9.5 


SS ODB 


0.5 


0.2 


9.6 


70.7 


5.6 


< 13.4 


ES ODB 


52.2 


14.1 


1.4 


1.6 


0.8 


< 29.9 




(55.9) 













TABLE II. Two-directional coupler: Efficiencies (in %) in 
addressing channels C22, C23, and Cn by an OVS, SS ODB, 
and ES ODB, respectively, and cross-talk efficiencies. 



Channel 
guiding by 


Cn 


C12 


Cl3 


C21 


C22 


C23 


C31 


C32 


C33 


Energy 
losses 


OVS 


1.3 


7.3 


1.3 


7.3 


61.6 


7.3 


1.3 


7.3 


1.3 


<4.0 


SS ODB 


0.2 


2.4 


12.6 


0.1 


6.1 


49.6 


0.2 


2.4 


12.6 


<13.7 


ES ODB 


47.1 


5.5 


0.7 


5.5 


0.3 


0.5 


0.7 


0.5 


0.5 


<38.6 
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